Optical excitation spectrum of an atom in a surface-induced potential 
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We study the optical excitation spectrum of an atom in the vicinity of a dielectric surface. We 
calculate the rates of the total scattering and the scattering into the evanescent modes. With a 
proper assessment of the limitations, we demonstrate the portability of the flat-surface results to 
an experimental situation with a nanofiber. The effect of the surface-induced potential on the 
excitation spectrum for free-to-bound transitions is shown to be weak. On the contrary, the effect 
for bound-to-bound transitions is significant leading to a large excitation linewidth, a substantial 
negative shift of the peak position, and a strong long tail on the negative side and a small short tail 
on the positive side of the field-atom frequency detuning. 

PACS numbers: 42.50.Vk,42.50.-p,32.80.-t,32.70. Jz 



The study of individual neutral atoms in the vicin- 
ity of dielectric and metal surfaces has gained renewed 
interest due to progress in atom optics and nanotech- 
nology 0, S H II II B 0. The ability of manipulating 
atoms near surfaces has established them as a tool for 
cutting-edge applications such as quantum computation 
H, 0, 13 , atom chips 0, 0], and probes, very sensi- 
tive to the surface-induced perturbations [Fj. Recently, 
translational levels of an atom in a surface-induced po- 
tential have been studied 

[111 US] An optical technique 
for loading atoms into quantum adsorption states of a 
dielectric surface has been suggested 0, Q . An experi- 
mental observation of the excitation spectrum of cesium 
atoms in quantum adsorption states of a nanofiber sur- 
face has been reported |3|. Spontaneous radiative decay 
of translational levels of an atom near a dielectric surface 
has been investigated Q. In this paper, we study the op- 
tical excitation spectrum of an atom in a surface-induced 
potential. 

We assume the whole space to be divided into two re- 
gions, namely, the half-space x < 0, occupied by a nondis- 
persive nonabsorbing dielectric medium (medium 1), and 
the half-space x > 0, occupied by vacuum (medium 2). 
We examine an atom, with an upper internal level e and 
a lower internal level g, moving in the empty half-space 
x > 0. 

The potential energy of the surface-atom interaction is 
a combination of a long-range van der Waals attraction 
— Cz/x 3 and a short-range repulsion Here C3 is the 
van der Waals coefficient. We approximate the short- 
range repulsion by an exponential function Ae~ ax , where 
the parameters A and a determine the height and range, 
respectively, of the repulsion. The combined potential 
depends on the internal state of the atom (see Fig. QJ, 
and is presented in the form Vj(x) — Aje~ ajX — C^j/x 3 , 
where j = e or g labels the internal state of the atom. 
The potential parameters Csj , Aj , and ctj depend on the 
dielectric and the atom. In numerical calculations, we 
use the parameters of fused silica, for the dielectric, and 
the parameters of atomic cesium with the D2 line, for 
the two-level atomic model. According to Ref. £|, the 



parameters of the ground- and excited-state potentials 
for the silica-cesium interaction are C 3g = 1.56 kHz /im 3 , 
C 3e = 3.09kHz ^m 3 , A g = 1.6xl0 18 Hz, A e = 3.17xl0 18 
Hz, and a g — a e = 53 nm -1 . 

We introduce the notation ip a = (p Ve and (ft, = (p v for 
the eigenfunctions of the center-of-mass motion of the 
atom in the potentials V e and V g , respectively. They are 
determined by the stationary Schrodinger equations 



2m dx 2 



+ V j (x) 



<p Vj (x) = £ Vs y Vj {x). 



(1) 



Here m is the mass of the atom. The eigenvalues £ a = £ Ve 
and £1, = £ v are the total center-of-mass energies of the 
translational levels of the excited and ground states, re- 
spectively. These eigenvalues are the shifts of the energies 
of the translational levels from the energies of the corre- 
sponding internal states. Without the loss of generality, 
we assume that the center-of-mass eigenfunctions ip a and 
(fib are real functions. 



FIG. 1: Energies and wave functions of the center-of-mass 
motion of the ground- and excited-state atoms in the surface- 
induced potentials. The parameters of the potentials are 
C 3g = 1.56 kHz /jm 3 , C 3e = 3.09 kHz /im 3 , A g = 1.6 x 10 18 
Hz, A e = 3.17 x 10 18 Hz, and a g = a e — 53 nm _1 . The mass 
of atomic cesium m = 132.9 a.u. is used. We plot, for the 
excited state, two bound levels (y — 400 and 415) and, for 
the ground state, two bound levels (y = 281 and 285) and 
one free level (£/ = 4.25 MHz). 

We introduce the combined eigenstates |o) = |e) (g> \ip a ) 
and I b) = \g) £g) \<fb), which are formed from the internal 
and translational eigenstates. The corresponding ener- 
gies are %uj a = htv e + £ a and Timb = hto g + £b- Here, ujj 
with j = e or g is the frequency of the internal level j. 
Then, the Hamiltonian of the atom moving in the surface- 
induced potential can be represented in the diagonal form 
Ha — Yla ^ w a|a)(a| + J2 b huJb\b)(b\. We emphasize that 
the summations over a and b include both the discrete 
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(£ a ,b < 0) and continuous (£ a ,b > 0) spectra. The levels 
a with £ a < and the levels b with £ b < are called the 
bound (or vibrational) levels of the excited and ground 
states, respectively. In such a state, the atom is bound 
to the surface. It is vibrating, or more exactly, moving 
back and forth between the walls formed by the van der 
Waals potential and the repulsion potential. The levels 
a with £ a > and the levels b with £ b > are called 
the free (or continuum) levels of the excited and ground 
states, respectively. The center-of-mass wave functions 
of the bound levels are normalized to unity. The center- 
of-mass wave functions of the free levels are normalized 
to the delta function of energy. For bound states, the 
center-of-mass wave functions ip a and ifb can be labeled 
by the quantum numbers v a and Vb, respectively. For 
free levels, the conventional sums over a and b must be 
replaced by the integrals over £ a and £b, respectively. 

Suppose that the atom is driven by a coherent plane- 
wave laser field Ei = [£ e~ t ^ kx+ut ^e + c.c.]/2, propagating 
perpendicularly to the surface of the dielectric. Here £ , e, 
u), and k = lu/c are the envelope, the polarization vector, 
the frequency, and the wave number, respectively, of the 
laser field. The field reflected from the surface is E2 = 
-R[£e i ( kx - Ut )e + c.c.]/2, where R = {n x - IV (ni + 1) is 
the reflection coefficient. According to Ref. |4j, the time 
evolution of the density matrix p of the atom is governed 
by the equations 

l,b 

- Tjila + 7a>)Paa>, 

P"* = \^2^lab'Pb'be~ l5ab,t ~ ^Y. nia ' b P aa ' e ~ lSa ' bt 
Lb' l,a' 
1 

— T^laPab-i 

PW = - l -Y.^pl b e- lS ^ t -^abPab'e^ t ) 

I. a 

+ \ X)(7a«W + Ja'ab'bV^'-^^Paa', (2) 
aa' 

where the parameters j a a>bb> and j a = J^b^abb are the 
radiative decay coefficients [4| , the index I = 1 or 2 labels 
the incident or reflected field, respectively, the notation 
5 a b = lo — Lu a + LUb stands for the detuning of the field 
frequency u> from the atomic transition frequency u) ab — 
uj a — ujb, and the quantities 

n lab = nF ab (-k) = QF* b (k), 

Q 2 ab = -RnF ab (k) (3) 

are the Rabi frequencies of the fields E;=i,2 with respect 
to the translational-state transition a <-» b. Here, — 
£d eg /h is the Rabi frequency of the driving field with 
respect to the transition between the internal states |e) 
and \g), with d eg = e ■ (e|d|g) being the projection of 



the atomic dipole moment onto the polarization vector 
e. The coefficient 

F ab (k) = (ip a \e ikx \<p b ) (4) 

is the overlapping matrix element that characterizes the 
transition between the center-of-mass states \ip a ) and \ip b ) 
with a transferred momentum of hk. Since (p a and (p b are 
real functions, we have F ab (—k) = F* b (k). In deriving 
Eqs. ©, we have neglected the cross decay coefficients 
7 aa ' = ^2 b ^aa'bb with a ^ a'. Such coefficients are small 
compared to the radiative linewidths j a and j a > when the 
refractive index n\ of the dielectric is not large 

We plot in Fig. [5]the Frank-Condon factors \F a b(k)\ 2 
for three bound ground-state levels v b and various bound 
excited-state levels v a . As seen from the figure, a deep 
bound ground-state level can be substantially coupled to 
a deep excited-state level. A deeper lower level can sub- 
stantially overlap with a deeper upper level. The level 
v b = 285, with translational energy shift £ b = —54.39 
MHz, can be substantially coupled to the level v a = 400, 
with translational energy shift £ a = —132.84 MHz. The 
difference £ a — £ b = —78.45 MHz is negative and large 
compared to the natural linewidth 70 = 5.25 MHz of 
the D2 line of atomic cesium. In principle, both types 
of shifts, namely red (negative) shifts (£ a — £ b < 0) and 
blue (positive) shifts (£ a — £b > 0), of the transition fre- 
quencies may be obtained. Note that the van der Waals 
potential for the excited state is stronger than that for the 
ground state. Therefore, for a bound-to-bound transition 
with a substantial Frank-Condon factor, the energy shift 
of the excited-state level is usually larger than that of 
the ground-state level. Consequently, the strong bound- 
to-bound transitions usually have red shifts. Since the 
depths of the surface-induced potentials V g and V e are 
large (159.6 THz for V g and 316 THz for V e ), the range of 
the frequency shifts of the strong bound-to-bound tran- 
sitions is broad, and the shifts can reach large negative 
values. 

We plot in Fig. [Slthe Frank-Condon factors \F a b f (k)\ 2 
for a free ground-state level b / and various bound excited- 
state levels a. The figure shows that the values of the 
Frank-Condon factors \F abf (k)\ 2 are more substantial for 
shallow bound excited-state levels than for deep ones. 
Among the bound excited-state levels, the level that is 
most strongly coupled to the free ground-state level 6/ 
is the level with the quantum number v a = 415. The 
translational energy shift of this level is £ a — —6.56 MHz. 
This bound level of the excited state is rather shallow. 
Due to this fact, the frequency shift £ a — £ bf — —10.81 
MHz of the strongest free-to-bound transition from the 
free level £ bf — 4.25 MHz is negative but not large. We 
note that, when the temperature of the atomic system 
is low, the typical values of the free-level energy £ bf are 
small. In this case, the range of the frequency shifts of 
substantial free-to-bound transitions is small compared 
to the range of the frequency shifts of substantial bound- 
to-bound transitions. 

We now calculate the excitation spectrum of the atom 
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FIG. 2: Frank-Condon factors |.Fa&(AOI for three bound 
ground-state levels b, with vb = 285, 290, and 295, and various 
bound excited-state levels a, with v a in the range from 385 to 
430. The light wavelength is A = 852 nm. Other parameters 
are as in Fig. 
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FIG. 3: Frank-Condon factors \F a b f (k)\ 2 for a free ground- 
state level bf, with £),, = 4.25 MHz, and various bound 
excited-state levels a, with v a in the range from 385 to 430. 
The parameters used for calculations are as in Fig. |2] 



in the surface-induced potential. Assume that the field is 
so weak that the Rabi frequency Q is much smaller than 
the natural linewidth 70. We also assume that the atom 
is initially in a coherent mixtures of the ground-state lev- 
els b with the weight factors pbb = Pb&(0), which are the 
initial populations of the levels. We consider the adia- 
batic regime. We find from Eqs. @ that the population 
p aa of the excited-state level a is given by 



E 



451 



\F* ab {k) - RF ab {k)\ 2 p bb , (5) 



where s = 2|f2| 2 /7g is the saturation parameter. The 
rate of total scattering of light from the atom is given by 



r sc = J2a Paaja- With the help of Eq. ©, we find 

2 

r - = i E ^ 1\*<2 l*\Kb(Q - RF ab (k)\ 2 Pbb . (6) 

ab ' a ab 

We note that, due to the presence of the interface, the 
field in a mode can be either a propagating light wave or 
a bound evanescent wave |l5l . An evanescent wave 
appears on the vacuum side of the interface when a light 
beam passing from the dielectric to the vacuum is totally 
internally reflected. Therefore, we can decompose the 
spontaneous emission rate as j a = 7^, + 7a j where 
ja ^ and ji 1 ^ are the rates of spontaneous emission 
from the level a into the bound evanescent modes and the 
propagating radiation modes, respectively 0. Hence, we 



have r sr = r 



(ev) 



^(ev) 



r 



it 



(rad) 



where 



2^ 

ab 



T^IKfcW-^WlW (7) 



is the rate of scattering from the atom into the evanescent 
modes and rsc ad ' ) = T sc — rsc V ' is the rate of scattering 
from the atom into the propagating radiation modes. 

All the above results were derived for a two-level atom 
moving in a potential induced by an infinite flat surface, 
and Eqs. @ and J7J represent the rigorous expressions 
for the rates of the total scattering and the scattering into 
the evanescent modes (in the framework of the adiabatic 
approximation and the perturbation theory). However, 
with a proper assessment of the limitations, the above 
theory can be applied to typical experimental situations. 
As an example, we pick a recent experiment , whereby, 
the number of photons scattered from cesium atoms into 
the guided modes of a nanofiber was measured as a func- 
tion of the frequency of the probe field. To apply our 
flat-surface results to the situation of a nanofiber, we use 
several assumptions and approximations. First of all, we 
note that, in the case of an atom near a fiber, the effective 
potential for the center-of-mass motion of the atom along 
the radial direction must include the centrifugal poten- 
tial U c { = Ti 2 {l1 — 1/4) /2mr 2 . Here l z is the quantum 
number for the axial component of the angular momen- 
tum of the atom. Since the van der Waals potential is 
t^vdw = C*3/r 3 , we have |J7 c f| <C |C/vdw| when the con- 
dition r < r c ee 2mC 3 /{h 2 \ll - 1/4 1) is satisfied. Thus, 
when the atom is near to the fiber surface and the axial 
angular momentum is small, the effect of the centrifugal 
is small compared to that of the van der Waals potential. 
In this case, we can neglect the centrifugal potential. For 
cesium atoms with m = 132.9 a.u. and l z = 10, the 
condition for the validity of the above approximation is 
r <C r c = 411 nm. Next, we note that the reflection of an 
incident plane wave from a cylindrical surface is rather 
complicated; it does not produce a plane wave. However, 
for the silica-vacuum interface, the reflection coefficient 
for light with wavelength A = 852 nm is R = 0.18, a 
small number. Therefore, to simplify our calculations, 
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we can neglect the reflection of light from the fiber sur- 
face. Finally, we note that there are two types of modes 
of the field in the presence of a nanofiber, namely guided 
modes and radiation modes. The guided modes are char- 
acterized by evanescent waves on the outside of the fiber. 
They are then the direct analog of the evanescent modes 
of a dielectric-vacuum interface. 

For a cesium atom at rest, the rate 7(r) of total sponta- 
neous emission and the rate 7^ (r) of spontaneous emis- 
sion into the guided modes of a nanofiber have been cal- 
culated systematically in Ref. @. For an atom in a 
translational state \a), the rate 7 a of total spontaneous 
emission and the rate 7^ of spontaneous emission into 
the guided modes can be estimated as Q 

/>oo 

7a = / j(r)\ip a (r)\ 2 dr, 
Jo 

poo 

7 fe) = / 7 te)(r)|^ a (r)| 2 dr. (8) 
Jo 

We plot in Fig. 21 the rates 7 a and 7 a s ' ) for a nanofiber 
with radius of 200 nm, which was used in the experiment 
. Figure 0] shows that the emission into the guided 
modes is substantial for a wide range of bound excited- 
state levels (v a < 430). However, when v a is very large 
iy a > 430), the channeling of emission into the guided 
modes becomes negligible. The reason is that, when v a 
is very large the bound state is very shallow. The atom 
in such state spends most of its time far away from the 
surface, where the evanescent waves in the guided modes 
cannot penetrate into. 
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FIG. 4: Rate of total spontaneous emission, j a (triangles), 
and rate of spontaneous emission into the guided modes, 7^ 
(quadrates) , for a nanofiber. The parameters for the surface- 
atom potentials are as in Fig. Q The atomic transition wave- 
length is Ao = 852 nm. The fiber radius is 200 nm. All the 
rates are averaged with respect to the orientation of the dipole 
moment. 

With the above simplifications and approximations, 
the rate of scattering into the guided modes of a nanofiber 



can be estimated by 

r^-fE^V^I^WPm. 0) 

The dependence of Tlf 1 on the field-atom detuning 5 = 
lu—luq (the difference between the probe field frequency oj 
and the free-space atomic resonance frequency loq) char- 
acterizes the optical excitation spectrum of the atom 
in the surface-induced potential. We calculate for 
two cases. In the first case, the atom is initially in free 
ground-state levels. In the second case, the atom is ini- 
tially in bound ground-state levels. We plot the results 
for the first and second cases in Figs. El and respec- 
tively. 

0.08 1 , 1 , , 1 1 
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FIG. 5: Rate Vsc of light scattering from the atom into the 
guided modes of a nanofiber through the free-to-bound tran- 
sitions as a function of the detuning S = lo — u>o of the probe 
field. The atom is initially in a thermal mixture of free ground 
states with temperature T = 200 fiK. The range of the quan- 
tum numbers v a of the bound excited states involved into the 
calculations is from 385 to 429. The natural linewidth of the 
atom is 70 = 5.25 MHz. Other parameters are as in Figs. 
andg| 

In the first case, i.e. the case of Fig. \5\ we assume that 
the initial state of the atom is a thermal mixture of free 
ground states \bf), described by the density operator 

/>oo 

p(0) = \g)(g\ P E \E)(E\dE. (10) 
Jo 

Here, \E) are the continuum eigenstates, normalized to 
the delta function of energy, for the ground-state Hamil- 
tonianp 2 / 2m + Vg(r)> and P E = e~ E/kBT jZ is the Boltz- 
mann weight factor, with T and Z being the temperature 
and partition function, respectively, for the atomic sys- 
tem. The dependence of the scattering rate on the 
field detuning 5 in Fig. ^represents the optical excitation 
spectrum for the atomic free-to-bound transitions. This 
spectrum shows a small negative shift of about —0.4 MHz 
for the position of the peak. The linewidth of the spec- 
trum is about 6.7 MHz, slightly larger than the natural 
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linewidth 70 = 5.25 MHz of the D2 line of atomic cesium. 
In addition, the spectrum has a small, short tail on the 
negative side of the detuning. These features are the re- 
sults of the negative frequency shifts of the transitions 
between free ground-state levels and bound excited-state 
levels. The observed effects are however weak, and the 
excitation spectrum is basically concentrated around the 
atomic resonance frequency loq . The reason is the follow- 
ing: The Frank-Condon factor for the overlap between a 
free ground-state and a deep bound excited-state level is 
small (see Fig. |3J). Therefore, a free ground-state atom 
can be excited only to shallow bound levels (and free lev- 
els) of the excited state. Meanwhile, an atom in a shallow 
(or free) excited level cannot emit photons efficiently into 
the guided modes because the atom spends most of its 
time far away from the fiber (see Fig. QJ. Hence, the 
effect of the surface-induced potential on the excitation 
spectrum in the case of free-to-bound transitions is rather 
weak. 
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FIG. 6: Rate Tif of light scattering from the atom into the 
guided modes of a nanofiber through the bound-to-bound 
transitions as a function of the detuning 6 = u) — (jJo of the 
probe field. The atom is initially in an incoherent mixture of 
a number of bound ground-state levels with quantum num- 
bers Vb from 269 to 293 and equal weight factors. The bound 
excited-state levels involved into the calculations are those 
with the quantum numbers v a from 385 to 429. The natural 
linewidth of the atom is 70 = 5.25 MHz. Other parameters 
are as in Figs. 0and^] 

In the second case, i.e. the case of Fig. we assume 
that the initial state of the atom is in an incoherent mix- 
ture of a number of bound ground-state levels Vb, with 
flat weight factors pbb = const. The initial density oper- 
ator of the atom is given by 

Vb= v xain 

where N is the number of bound ground-state levels in- 
volved in the initial state of the atom. The use of flat 
weight factors is explained by the assumption that the 



atom is adsorbed by a surface at the room temperature 
T = 300 K = 6.25 THz, which is much higher than the 
energy differences between the levels. For our numerical 
calculations, we include the bound ground-state levels 
with quantum numbers Vb from 269 to 293. The energy 
shifts of these levels, equal to their total center-of-mass 
energies, are in the range from —I GHz to —5 MHz. We 
do not include levels with larger quantum numbers be- 
cause they are too shallow and therefore, an atom in such 
a bound ground-state level can be easily excited to a free 
ground-state level by heating or collision. We also do not 
include levels with smaller quantum numbers because the 
corresponding transition frequency shifts are beyond the 
range of the excitation spectrum measured in the exper- 
iment 0- The dependence of the scattering rate vi^ 
on the field detuning 6 in Fig. [B] represents the opti- 
cal excitation spectrum for the atomic bound-to-bound 
transitions. This spectrum shows a substantial negative 
shift of about —14.3 MHz for the position of the peak. 
The linewidth of the spectrum is about 58.3 MHz, which 
is one order larger than the natural linewidth 70 = 5.25 
MHz of the D2 line of atomic cesium. Furthermore, the 
spectrum has a substantial long tail on the negative side 
and a small short tail on the positive side of the detun- 
ing. These features are the results of the frequency shifts 
of the transitions between bound ground-state levels and 
bound excited-state levels. In general, both negative and 
positive shifts can come into play. However, the transi- 
tions with positive shifts are relatively weak because the 
corresponding Frank-Condon factors are usually small. 
The overlap between the center-of-mass wave functions 
of a bound excited-state level and a bound ground-state 
level is substantial only when the corresponding return- 
ing points are close to each other. Since the van der 
Waals potential for the excited state |e) is deeper than 
that for the ground state \g), the frequency shifts of the 
transitions that are associated with substantial Frank- 
Condon factors are mostly negative. The deeper the 
bound ground-state levels involved in the initial state, 
the stronger the effect of the surface-atom interaction 
on the excitation spectrum. We note that the basic fea- 
tures of the spectrum in Fig. [5] are very similar to the 
features of the right and left wings of the experimental 
spectrum reported in Ref. ,3]. However, it is impor- 
tant to recall that the experimental spectrum Q has con- 
tributions from both bound-to-bound and free-to-bound 
types of transitions. A combination of both Figs. [5] and 
[fj] with proper weight factors can indeed lead to an ex- 
cellent matching of the experimental spectrum with the 
theoretically predicted one. 

In conclusion, we have studied the optical excitation 
spectrum of an atom in the vicinity of a dielectric surface. 
We have derived rigorous expressions for the rates of the 
total scattering and the scattering into the evanescent 
modes. With a proper assessment of the limitations, our 
theoretical results have been applied to an experimental 
situation with a nanofiber. We have shown that the effect 
of the surface-induced potential on the excitation spec- 
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trum in the case of free-to-bound transitions is rather 
weak. Meanwhile, the spectrum of excitation of bound- 
to-bound transitions has a large linewidth, a substantial 
negative shift of the peak position, and a substantial long 
tail on the negative side and a small short tail on the posi- 
tive side of the field-atom frequency detuning. The basic 
features of the calculated spectra are in agreement with 
the experimental observations 
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